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ABSTRACT 


For  rapid  crack  propagation  in  an  elastic  perfectly-plastic  material, 
explicit  expressions  have  been  obtained  for  the  dynamic  strains  on  the 
crack  line,  from  the  moving  crack  tip  to  the  moving  elastic-plastic 
boundary.  The  method  of  solution  uses  power  series  in  the  distance  to 
the  crack  line,  with  coefficients  which  depend  on  the  distance  to  the 
crack  tip.  Substitution  of  the  expansions  in  the  equations  of  motion, 
the  yield  condition  (Huber-Mises)  and  the  stress-strain  relations,  yields 
a  system  of  nonlinear  ordinary  differential  equations  for  the  coefficients. 
These  equations  are  exactly  solvable  for  Mode-Ill,  and  they  have  been 
solved  in  an  approximate  manner  for  Mode-I  plane  stress.  The  crack-line 
fields  have  been  matched  to  appropriate  elastic  fields  at  the  elastic- 
plastic  boundary.  For  both  Mode-Ill  and  Mode-I  plane  stress,  the  plastic 
strains,  which  depend  on  the  elastodynamic  stress  intensity  factor  and  the 
crack-tip  speed,  have  been  used  in  conjunction  with  the  crack  growth 
criterion  of  critical  plastic  strain,  to  determine  the  relation  between 
the  far-field  stress  level  and  the  crack-tip  speed. 
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1.  Introduction 


At  high  crack-tip  speeds  the  mass  density  of  a  material  affects  the 
fields  of  stress  and  deformation  in  the  vicinity  of  a  propagating  crack 
tip.  For  essentially  brittle  fracture,  near-tip  dynamic  effects  have 
been  investigated  extensively  on  the  basis  of  linear  elastic  fracture 
mechanics.  By  now,  several  papers  have  reviewed  the  computation  of 
elastodynamic  stress  intensity  factors,  and  they  have  discussed  dynamic 
effects  on  the  fracture  criterion  of  the  balance  of  rates  of  energies, 
see  Achenbach  (1) ,  Freund  (2)  and  Kanninen  (3) .  The  combined  effect  of 
plastic  deformation  and  mass  density  on  near-tip  fields  has  not  yet 
received  that  much  attention.  This  is  not  surprising,  considering  the 
difficulties  that  are  encountered  in  the  quasi-static  analysis  of  fields 
near  a  growing  crack  in  an  elastic-plastic  material. 

For  quasi-statically  growing  cracks  the  asymptotic  structure  of 
near-tip  fields  in  elastic  perfectly-plastic  solids  has  been  analyzed  in 
considerable  detail.  A  recent  review  by  Rice  (4)  includes  a  general 
formulation,  and  it  presents  detailed  results  for  isotropic  materials  of  the  Huber- 
Mises  type.  In  general,  the  analytical  near-tip  results  must  be  supple¬ 
mented  by  numerical  calculations  to  determine  certain  arbitrary  functions 
that  appear  in  the  asymptotically  valid  near-tip  results. 

In  recent  papers,  Achenbach  and  Dunayevsky  (5)  and  Achenbach  and  Li  (6) 
have  constructed  quasi-static  solutions  that  are  valid  on  the  crack  line, 
from  the  moving  crack  tip  up  to  the  moving  elastic-plastic  boundary. 
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These  solutions  were  obtained  for  an  elastic  perfectly-plastic  material  of 
the  Huber-Mises  type  by  expanding  all  fields  in  powers  of  the  distance, 
y,  to  the  crack  line.  Substitution  of  the  expansions  in  the  equilibrium 
equations,  the  yield  condition  and  the  constitutive  equations  yields  a 
system  of  simple  ordinary  differential  equations  for  the  coefficients 
of  the  expansions.  As  shown  in  (6),  the  resulting  equations  are  exactly 
solvable  for  the  Mode-Ill  case,  and  they  are  solvable  for  the  Mode-I  plane- 
stress  case  if  it  is  assumed  that  the  cleavage  stress  is  uniform  on  the 
crack  line.  By  matching  the  relevant  stress  components  and  particle 
velocities  to  the  dominant  terms  of  appropriate  elastic  fields  at  the 
elastic-plastic  boundary,  the  plastic  strains  on  the  crack  line  were 
computed  in  terms  of  the  elastic  stress  intensity  factor. 

The  literature  on  dynamic  effects  in  the  presence  of  elastic-plastic 
constitutive  behavior  is  growing.  Investigations  of  the  asymptotic 
structure  of  the  dynamic  near-tip  fields  were  presented  by  Slepyan  (7) 
and  Achenbach  and  Dunayevsky  (8).  Dynamic  near-tip  effects  for  a  strain¬ 
hardening  material  were  investigated  by  Achenbach  and  Kanninen  (9)  and 
Achenbach,  Kanninen  and  Popelar  (10)  on  the  basis  of  J^-flow  theory  and 
a  bilinear  effective  stress-strain  relation.  For  Mode-Ill  crack  propa¬ 
gation  in  an  elastic  perfectly-plastic  material,  exact  crack-line  solutions 
were  obtained  by  Achenbach  and  Dunayevsky  (11)  and  Freund  and  Douglas  (12) . 

In  the  present  paper  the  expansion  technique  of  Achenbach  and  Li  (6) 
is  extended  to  the  dynamic  formulation,  for  rapid  crack  growth  in  Mode-Ill 
and  in  Mode-I  plane  stress.  Systems  of  nonlinear  ordinary  differential 
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equations  have  been  established  which  are  valid  for  the  transient  case. 
Solutions  have,  however,  been  obtained  only  for  the  steady-state  dynamic 
crack  line  fields.  The  equations  for  the  Mode-Ill  case  can  be  solved 
rigorously  in  implicit  form.  An  approximate  approach  which  gives 
excellent  results  for  the  Mode-Ill  case  has,  however,  also  been  developed. 
The  equations  for  Mode-I  plane  stress  cannot  be  solved  rigorously,  but 
the  approximate  approach  can  be  used  to  yield  the  steady-state  dynamic 
cleavage  strain  on  the  crack  line.  The  plastic  strains  on  the  crack  line 
have  been  used  in  conjunction  with  the  crack  growth  criterion  of  critical 
plastic  strain  to  determine  a  relation  between  the  far-field  stress 
level  and  the  crack-tip  speed. 

The  geometry  is  shown  in  Fig.  1.  The  x^-axis  of  a  stationary  coordinate 
system  is  parallel  to  the  crack  front,  and  x^  points  in  the  direction  of 
crack  growth.  The  position  of  the  crack  tip  is  defined  by  x^  *  a(t). 

A  moving  coordinate  system  (x,y,z)  is  centered  at  the  crack  tip,  with  its 
axes  parallel  to  the  x^*x2  and  x^  axes. 
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2.  Mode  III  Crack  Propagation 

In  this  Section  an  exact  steady-state  dynamic  solution  is  derived 
which  is  valid  on  the  crack  line  in  the  plastic  loading  zone  ahead 
of  the  propagating  crack  tip. 

In  the  moving  coordinate  system  the  equation  of  motion  is 


3x  3x 
. xz  _jrz 

3x  3y 


pv 


(2.1) 


where  w(x,y,t)  is  the  anti-plane  displacement,  and  the  material  time 
derivative  is  defined  as 

C)  -  (3/3t)  -  a  (3/3x) 

Here  a  *  da/dt  is  the  speed  of  the  crack  tip.  The  Huber-Mises  yield 

condition  requires 

t2  +  t2  ■  k2  , 
xz  yz 

where  k  is  the  yield  stress  in  pure  shear.  The  strain  rates  are 


(2.2) 


(2.3) 


1  3w 


1  3w 


xz  2  3x  ’  yz  2  3y 

The  strain  rates  are  related  to  the  stresses  and  the  stress  rates  by 

I  -Ia+  At  ,  i  at 

xz  2  p  xz  yz  2  p  yz 

In  (2.5a,b)  p  is  the  shear  modulus  and  A  is  a  positive  function  of 
time  and  the  spatial  coordinates. 

Solution  along  the  crack  line. 

In  this  paper  we  are  interested  in  solutions  along  the  crack  line 
y  ■  0,  0  <  x  Xp,  where  x  ■  x^  defines  the  elastic-plastic  boundary. 
Such  solutions  can  be  obtained  by  considering  expansions  with  respect 
to  y  in  the  region  y/x  <<  1: 


(2.4a,b) 


(2.5a,b) 


4 


(2.6) 


i  =  so(x,t)  +  s,,(x,t)y2  +  0(y4) 
Txz  =  T1(x,t)y  +  0(y3) 
w  =  w1(x,t)y  +  0(y3) 


A  =  A  (x,t)  +  A  (x,t)y  +  0(y4) 

O  2 

In  (2. 6) -(2. 9)  we  have  taken  into  account  that  t  and  A  are 
symmetric  with  respect  to  y  ■  0,  while  t  and  w  are  antisymmetric. 
Substitution  of  (2.6)-(2.9)  into  (2.1),  (2.3)  and  (2.5a,b),  and 
collecting  terms  of  the  lowest  orders  in  y  yields 
3t, 


lt+2s2 


P«, 


s2  *  k2  ,  2s  s-  +  t2  =  0 
o  o  2  i. 


1  3*1  fi  1.  ' 

—  — r—  *  —  +  At,,  -=w,  *=  A  s 

2  8x  2  m  ol  21  oo 


It  follows  from  (2.11a)  that  s  ■  k.  Elimination  of  s„  from  (2.10) 

o  2 

and  (2.11b)  gives 
3t  tJ 

-si  -  r  -  p*i  - 0 

Similarly,  Aq  can  be  eliminated  from  (2.12a)  and  (2.12b)  to  yield 

!l  I  •  ,  n 

3x  "  ,  ’  k  11  "  0 

Equations  (2.13)  and  (2.14)  define  two  coupled  nonlinear  partial 
differential  equations.  Analytical  solutions  to  these  equations, 
which  would  give  the  transient  fields  on  the  crack  line,  have  not 
yet  been  obtained. 


(2.7) 

(2.8) 

(2.9) 

(2.10) 

(2.11a»b) 

(2 ,12a, b) 

(2.13) 

(2.14) 
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Equations  (2.13)  and  (2.14)  must  be  supplemented  by  conditions 

at  the  elastic-plastic  boundary  E.  These  have  been  discussed  in  some 

detail  in  Appendix  A,  where  it  was  shown  that  for  conditions  which  may 

be  assumed  to  hold  ahead  of  a  propagating  crack  tip,  the  stresses  are 

continuous  at  E,  see  Eq. (A.24a,b) .  From  the  impulse  momentum  relation 

(A. 3)  it  follows  that  the  particle  velocity  is  then  also  continuous. 

Near  the  crack  line  at  x  ■  x  we  can  then  write: 

P 

t  [sq]  ]  -  0  ,  [[t1]]  -  0  and  [[w^]  =  0  ,  (2.15a,b, 

where  the  notation  for  discontinuities  is  defined  by  Eq.(A.l). 

The  governing  equations  for  the  quasi-static  case  follow  by 
setting  p  =  0.  The  resulting  system  of  coupled  nonlinear  ordinary 
differential  equations  can  be  solved.  The  quasi-static  solution  for 
w^  has  been  given  in  Ref. (5). 

For  the  steady-state  case  the  material  time  derivative  (2.2) 
reduces  to 

C)  -  -  a  (d/dx)  (2.16) 

where  a  is  now  a  constant  crack  tip  speed.  We  define 


where  the  Mach  number  M  is  defined  as 


M  =  a/(u/p)'1. 


(2.20) 


For  small  values  of  y  (i.e.,  y/x  <<  1)  the  plastic  fields  in  the 
loading  zone  will  be  matched  to  the  dominant  terms  of  the  elastic 
fields.  In  polar  coordinates  R,i|/  centered  at  point  E,  and  for  small 
values  of  the  angle  \p,  the  dominant  terms  of  the  solution  on  the  elastic 
side  of  the  elastic-plastic  boundary  are  taken  as 


Here  the  elastic  stress-intensity  factor  depends  on  M.  The 

angular  dependence  on  M  enters  in  higher  order  terms  of  <|/.  It  should 
be  noted  that  the  center  of  the  elastic  field  is  not  taken  to  coincide 
with  the  crack  tip.  The  center  is  located  at  a  moving  point  E.  The 
geometry  is  shown  in  Fig.  1. 

Since  is  continuous  at  y  *  0,  x  *  xp,  see  Eq. (2.15a),  we  find 
(lirir)  kiii  “  k  ’  or  RP  *  (Km>2/27rk2 


(2.21) 

(2.22a,b) 


(2.23a,b) 


where  R  =  R^  defines  the  radius  of  curvature  of  the  elastic-plastic 
boundary,  at  least  for  small  values  of  \p.  Another  condition  is  that 
(i.e.,  the  shear  stress  in  the  R,ij>  system)  should  be  continuous 
at  the  elastic-plastic  boundary.  We  find  by  using  (2.6)  and  (2.7) 


tRZ  “  Txz  C08,,'  +  Tyz  8in<'  S  Ti7  +  k* 


(2.24) 
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Thus,  by  the  use  of  (2.22a)  and  (2.23a) 


Tjy  +  k<J>  = 

Since  ip  s  y/R  ,  we  obtain  at  x  =  x 
P  P 


t  =  -  k/2R 
1  P 

Continuity  of  y.  at  x  *  x  ,  which  follows  from  (2.15c),  yields 
1  P 


3w  3w  1,, 
- x  —  =  -  -^Jjk 


Yly  3R  R  dip 


or 


1 


(2.25) 


(2.26) 


(2.27) 


(2.28) 


Y,  =  ~  k/2yR 
1  P 

where  (2.23a)  has  been  used.  For  completeness  we  list  the  condition 
on  the  strain  3w/3y  at  x  *  x^ 

Y  =  I2  =  « .  =  -  (2.29) 

y  3y  1  y  v  ' 

Equations  (2.18)  and  (2.19)  can  be  solved  rigorously,  as  shown 

in  Appendix  B.  It  is,  however,  of  interest  to  note  that  an  asymptotic 

solution  for  small  values  of  x  can  be  obtained  with  minimal  effort. 

Let  us  consider  solutions  of  the  general  form 


^  3  -  Tj/x  , 


Y1  »  V7* 


(2.30a,b) 


Substitution  in  (2.18)  and  (2.19)  yields 
-  k(l±M)  ,  Y  x  m  -  (im)k/VM 

Since  we  must  have  y^  <  0,  we  discard  the  solution  containing  the 
plus  signs.  Hence 


(2.31a,b) 


8 


k(l-M) 


k  I;4t  l 
V  M  x 


.  32a,b) 


Since  Yj  =  9Yy/3x,  we  also  have 
_  k  1-M  .  .  ,  . 

3  “  ~  -TT  *.n(x/x  )  (/ 

y  y  M  p 

This  solution  is  the  same  as  the  one  derived  earlier  by  Slepyan  (7), 

see  also  Achenbach  and  Dunayevsky  (8) .  Note  that  reduces  to  the 

quasi-static  solution  as  M  -*■  0.  The  strain  has,  however,  not  only 

the  wrong  behavior  in  x,  but  actually  becomes  singular  in  M. 

As  shown  in  Appendix  B,the  solution  to  Eqs.(2.18)  and  (2.19)  which 

satisfies  the  boundary  conditions  (2.26)  and  (2.28)  at  x  =  x^  is 

defined  by  the  following  equations: 

/  ,  vO-M)/2M 


/  1  V  ' 1 

T1  /1-M\1/2M  M  (~2F  ~  1+M/ _ 

k/2R  -“UW 

' 2F  +  I^m/ 


where 


b-  dt 


T  -I  ,  dT 

-  —  (-k(l-M2)  —  ~  +  1] 
UM2  dx 


Equation  (2.34)  gives  F  as  a  function  of  tj.  Integration  of  (2.35) 
then  yields 

-k/2R  52F(0  P 
P 
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Equation  (2.37)  yields  as  a  function  of  x  and  x^.  Substitution  of 

the  result  in  Eq.(2.36)  yields  y^.  By  letting  -*•  *  in  (2.37)  we  obtain 

a  relation  between  x  and  R  .  From  y,  =  3y  /i)x,  we  finally  find 

P  P  1  y 

x 

y  =  -  +  /  y,  dx  =  ±  +  yP  (2.38) 

y  v  Jx  1  w  y 

p 

The  strain  y^  obtained  from  (2.38)  is  the  exact  solution  on  the  crack 
line.  This  solution  is  equivalent  to  the  one  obtained  earlier  by 
Dunayevsky  and  Achenbach  (8)  ,  and  Freund  and  Douglas  (12).  It  can 
be  shown  that  for  small  x,  Eq.(2.38)  reduces  to  (2.31b).  In  the 
limit  M  -*  0,  (2.38)  reduces  to  the  quasi-static  solution 

~y  -  1  -  S.n(x/x  )  +  ~  («.n(x/x  )]2  (2.39) 

k  y  p  z  P 

An  explicit  analytical  expression  for  y^,  albeit  an  approximate 
one,  would  be  very  useful  for  applications  in  conjunction  with  the 
crack  growth  criterion  of  a  critical  plastic  strain.  Another  reason 
for  an  approximate  approach  to  the  Mode-Ill  case  is  that  the  results 
can  be  tested  by  comparison  with  exact  results.  The  same  approach 
can  then  be  used  for  the  Mode-I  plane-stress  case,  which  is  not 
amenable  to  an  exact  solution. 

An  approximate  approach  is  suggested  by  the  structure  of  Eqs.(2.18) 
and  (2.19).  If  an  acceptable  approximation  to  would  be  available 
a-priori,  then  (2.18)  would  simply  be  a  linear  ordinary  differential 
equation  for  y^.  A  first  approximation  to  is  suggested  by  (2.32a), 
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namely,  =  -k(l-M)/x.  This  expression  has  the  correct  limits  at 
M  =  0  and  M  =  1.  A  better  result  is  obtained  by  adding  a  constant 
term 

Tj_  -  -  k(l-M)[i  +  —-]  (2.40) 

P 

The  second  term  is  chosen  so  that  (2.19)  is  satisfied  up  to  order 
0(M)  near  x  =  xp.  Figure  2  shows  a  comparison  between  (2.40)  and 
the  exact  result.  Equation  (2.40)  can  now  be  substituted  in  (2.18), 
and  the  resulting  equation  can  be  solved  rigorously  for  y^.  The 
strain  y^  =  k/p  +  y^  then  follows  from  (2.38).  In  anticipation  of 
difficulties  with  the  Mode-I  case,  we  elect,  however,  to  solve  y^  by 
using  a  perturbation  solution  which  ignores  terms  of  order  0(M2). 

The  corresponding  expression  for  yy  is  obtained  as 

PY  P,x,  (1-M)  (2-3M-M2)  [l,x_.M  1  M(l-M)  .x  .1'H<  M(l-M) 
k  ^  Xp  "  M(2-M+M2)  [M  Xp  "  M  "  2(14M>  xp  MHO' 

-  ~  +  i(l-M)(^-  -  1)  (2.41) 

P  P 

It  is  of  interest  that  (2.41)  yields  (2.33)  in  the  limit  x  -*■  0, 
while  it  yields  the  quasi-static  solution  (2.39)  as  M  ■+  0, 
provided  that 

M  *n(^-)  <<  1  (2.42) 

XP 

A  comparison  of  (2.41)  with  the  exact  result  is  shown  in  Fig.  3. 
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Application  of  the  boundary  condition  (2.26)  to  the  approximate 


expression  for  t^,  yields  a  relation  between  xp  and  Rp .  Subsequent 
use  of  (2.23b)  gives 

xp  =  (1-M)  (1+2M)  ~  (Km/k)2 

In  Ref. (12),  Freund  and  Douglas  use  numerical  results  of  a  finite 

element  analysis  to  derive 

x  =  (0.295  -  0.5M2)(K  T/k)2 

p  Ill 

For  arbitrary  M  we  cannot  obtain  an  exact  analytical  expression  for  xp 
as  a  function  of  Kj^/k  from  (2.34)-(2.37)  .  However,  as  M  +  0,  we  find 

%  ‘  <2-6"2)  K<Km'k,! 

Figure  4  shows  a  comparison  between  (2.43),  (2.44)  and  the  exact 

relation,  which  follows  from  (2.34)-(2.37) . 

Finally,  following  Freund  and  Douglas  (12)  we  apply  the  crack  growth 

criterion  of  critical  plastic  strain  to  determine  the  value  of  Kjjj  that 

would  be  required  for  crack  growth  at  a  given  value  of  M.  The  crack- 

growth  criterion,  originally  proposed  by  McClintock  and  Irwin  (13), 

states  that  the  crack  will  grow  with  (normalized)  plastic  strain 

(u/k)Yy  =  Yj  at  x  *  Xj  on  y  *  0.  For  plastic  strain  below  at  x  *  x^ 

the  crack  cannot  grow.  As  discussed  by  Rice  (14)  the  characteristic 

length  x  is  related  to  K  ,  the  value  of  the  Mode  III  stress  intensity 
i  c 

factor  which  is  required  to  satisfy  the  fracture  criterion  for  a  stationar; 
crack,  by  the  relation 


(2.43) 


(2.44) 


(2;45) 
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(2.46a  ,b] 


I 


I 


3.  Mode-I  Crack  Propagation  in  Plane  Stress 


We  consider  a  state  of  generalized  plane  stress,  hence  a^,  o ^  and 

o  vanish  identically.  Relative  to  the  moving  coordinate  system  the 
yz 

equations  of  motion  are 


3a  3t 
_J£  +  „ 

3x  3y 


pu 


3t  3a 

+  _J 

3x  3y 


pv 


The  Huber-Mises  yield  condition  becomes 


a2  +  a2  -  a  a  +  3t2 


3k2, 


y  x  y  xy 
where  k  is  as  in  Eq.(2.3).  The  strain  rates  are 


(3.1a,b) 


(3.2) 


3u  .  m  dv  .  _  l/3u  3v\ 
3x  *  Gy  3  3y  *  Cxy  =  2\3y  3x/ 


y  3y  xy  2\3y 
The  strain  rates  are  related  to  the  stresses  and  stress  rates  by 


=  4(°  ~  vd  )  ■*•  4  A(2 a  -  a  ) 

3x  E  x  y  3  '  x  y 


~  -  4(°  -  vo  )  +  4  A(2o  -  o  ) 

3y  E  y  x  3  y  x 

l/3u  .  3v\  1+v  .  ’ 

2\3y  3x/  E  xy  xy 

where  E  and  v  are  Young's  modulus  and  Poisson's  ratio,  respectively, 
and  A  is  a  positive  function  of  time  and  the  spatial  coordinates. 

Solution  along  the  crack  line 

Analogously  to  (2.6)-(2.9)  we  consider 

ox  -  PD(x,t)  +  P2<x»t;)y2  +  p4(x,t)y1'  +  0(y6) 


(3.3a,b,c) 


(3.4) 

(3.5) 

(3.6) 


(3.7) 


14 


(3.8) 


°y  “  +  q2^x,t^y2  +  Vx,t)y4  +  °<y6> 

T  =  s1(x,t)y  +  s3(x,t)y3  +  0(y5) 


u  *  uQ(x,t)  +  u2(x,t)y2  +  0 (y1* ) 
v  *  v1(x,t)y+  v3(x,t)y3  +  0(y5) 


A  *  A  (x »t)  +  A  (x,t)y2  +  0(y4) 
o  l 

Here  we  have  taken  into  account  that  o  ,  o  ,  u  and  A  are  symmetric  with 

x  y 

respect  to  y  =  0,  while  t  and  v  are  antisymmetric.  Substitution  of 
(3.7)-(3.9)  into  (3.1a,b)  and  collecting  terms  of  the  same  order  in  y 
yields 


o  , 

- —  +  s.  =  pu 
3x  1  o 


3P2 

ST  +  3S3  =  PU2 


3sl  ..  3s3 

3^r  +  2q2  “  - 


PV1  *  3T  +  4q4  "  PV3 


Substitution  of  (3.7)-(3.9)  into  the  yield  condition  (3.2)  yields  by  the 
same  procedure 


p2  +  q 2  -  p  q  =■  3k2 
o  o  o  o 


<2po'qo)p2  +  (2VPo)q2  +  3sl  "  0 

p2  +  <2po’qo>P4  +  q2  +  (2VPo)q4  "  P2q2  +  6sl83  "  0 
Another  5  equations  are  obtained  by  using  (2.16)  and  (3.7)-(3.10)  in 
(3.4)-(3.6).  These  equations  have  been  listed  as  Eqs. (21)-(25)  by 


(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13a,b) 

(3.14a,b) 

(3.15) 

(3.16) 

(3.17) 
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Achenbach  and  Li  (6) ,  and  they  are  not  reproduced  here. 

At  this  stage  we  have  14  unknowns  and  12  equations.  Clearly,  the 
system  cannot  be  solved  without  further  simplifying  assumptions.  For  the 
quasi-static  problem  (i.e.,  p  5  0),one  assumption,  namely  that  qQ  *  constant, 
suffices  to  produce  a  solvable  system  of  equations,  as  shown  by  Achenbach 
and  Li  (6).  For  reference  purposes  we  state  the  coefficients  for  the 
quasi-static  stresses  obtained  in  (6) 


k  ,  P2 
2k  ,  q. 


3  k_ 
2x2 


s  =  0  ,  s  - - 

1  3  3 

X 

An  approximate  solution  can  be  obtained  for  the  steady-state  dynamic 
problem,  when  (’)  *  -ad/dx.  Equations  (3.13)-(3.14)  then  become 


(3 

(3 

(3 


dpo  . 
dx  1 


ds. 


EM 


d2u 
_ o 

dx2 

d2v, 


dp2 

~  +  3s_ 
dx  3 


EM 


d2u 

2  _ 2 


dx" 


dx"  +  2q2  * 


dx" 


ds. 

dx~ 


+  4q, 


EM 


d2v„ 


dx2 


where  M  is  defined  as 
M  -  a/ (E/p)*5 

The  equations  that  can  be  obtained  from  (3.4)-(3.6)  now  become 
d2u 
dx2 


!uo  UdVo  dqo\^l  dA0„ 

Z  ”  E\dx  "  V  dx  /  +  3  dx  (2po“qo^ 


(3 


(3 


(3 


(3 


. 18a,b) 

.19a,b) 
. 20a,b) 

.21a, b) 

.22a,b) 

.23) 

.24) 
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1/^2  wl\, 

E  V  dx  ~  dx  )+  3  dx  (2p2“q2)  +  3  dx  ^o^c 


(3.25) 


(3.26) 


(3.27) 


(3.28) 


^1  _  lA>  _  .  . 

dx  E\dx  V  dx  )  3  dx  qo  Po 

dvi  1  /dq?  dP,\  1  dA„  1  dA? 

3  dx~  =  E  (dx""  "  V  dx“)  +  I  dT"(2q2  "  P2)  +  3  dhT(2qo  '  Po5 

du.  .  d2v,  ,  .  ds,  dA 

_ 2  1  _ 1  _  1+v  1  _ o 

dx  2  2  E  dx  dx  1 

dx 

To  determine  a  solution  to  Eqs . (3.15)— (3.17) ,  (3 ,21)-(3 .22)  and 
(3.24)-(3.28) ,  we  start  by  making  the  same  assumption  as  for  the  quasi¬ 
static  case,  namely,  that  qQ  ■  constant.  A  second  assumption  is  that 
2pQ  -  qQ  *  c,  where  e  *  e(M),  but  e  <<  k.  It  then  follows  from  (3.15) 


PQ  -  k  +  0(e),  qQ  *  2k  +  0(e2),  2qQ-po  =  3k  +  0(e),  2pQ-qo  =  0(e)  (3 .29a ,b,c ,d) 

Since  both  pQ  and  qQ  are  constant,  Eq.(3.26)  implies  that  dAQ/dx  * 

(l/k)dv1/dx  +  0(e),  and  it  subsequently  follows  from  (3.24)  that 

d2u  /dx2  »■  0(e).  Substitution  of  these  results  in  (3.21a)  gives  s,  =  0(eM2), 

while  (3.22a)  gives  q ^  »  0(M2).  Next,  we  conclude  from  (3.16)  that 

+  3kq^  +  O^M1*)  -  0,  which  implies  that  e  *  0(M2).  Application  of  the 
preceding  results  to  (3.17)  gives  q^  ■  -(l/3k)p2  +  0(M2).  Substitution  of 

the  latter  result  in  (3.22b),  and  then  in  (3.21b)  yields 


24?  ?  u2  ? 

“•  +  £  ?2  "  EM — —  +  0(M2 ) 


(3.30) 
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Note  that  the  inertia  term  has  not  been  neglected  in  this  equation. 


since  it  provides  a  coupling  with  equations  for  u^* 
s1  -  0(eM2)  -  0(M4)  into  (3.28)  yields 


ooo 


Substitution  of 


Finally,  by  using  (3.31),  as  well  as  (3.29)  and  dA  /dx 

o 

Eq.(3.25)  gives 


2^r+  E3T+3kdT  P2  *  0(M  > 


(l/k)dv1/dx. 


(3.31) 


(3.32) 


In  a  further  reduction  we  ignore  the  terms  of  0(M2),  and  we  eliminate 
u2  by  the  use  of  (3.31),  to  obtain 


£fl+*-2  + 

1  TTW  2 

d3v 

lx 

—  EM 

■  0 

dx3 

(3.33) 

1  .  2 

4.  1  dP2  n 

2  dx2  3k  F 

*2  Vlx 

E  dx  *  0  ’ 

(3.34) 

where 

vlx  "  dvi/dx  (3.35) 

Equations  (3. 33) -(3. 34)  will  be  used  to  analyze  the  Mode-I  plane-stress 
fields. 

The  solutions  for  vlx  and  p2  oust  satisfy  certain  conditions  at  the 
elastic  plastic  boundary  E.  In  Appendix  A  it  was  shown  that  for  conditions 
which  may  be  assumed  to  hold  ahead  of  a  propagating  crack  tip ,  the 
stresses  are  continuous  at  E  ,  see  Eq.(A.31).  From  the  impulse  momentum 
relation  (A. 3)  it  then  follows  that  the  particle  velocity  is  also  continuous 
at  E. 
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For  y/x  <<  1  the  plastic  fields  in  the  loading  zone  will  now  be 
matched  to  the  dominant  terms  of  the  elastic  field.  In  polar  coordinates 
centered  at  point  E,  the  elastic  field  for  small  values  of  iji  is  taken 
as 


a  = 
x 


^2irR^  KI{(1  '  2R  (8  ~  16R^2} 


(3.36) 


1 


9n 


^2irR^  Ki^1  +  2RJ  +  (8  16R^  } 


(3.37) 


xy 


4ttR;  i'2  4R'^ 


(3.38) 


-  1  +  s>  +  i<5  -« ->2> 


(3.39) 


^  k  Ki<K  - 1  +  ■ 


(3.40) 


where  <  *  (3-v) / (1+v) .  This  elastic  field  has  one  more  parameter,  namely  n. 


than  the  usual  elastic  crack-tip  field.  Equations  (3.36)  -  (3.40) 


actually  correspond  to  the  field  for  a  notch  with  -jn  as  its  tip-radius 


of  curvature,  see  Creager  and  Paris  (15).  The  elastic  stress-intensity 
factor  Kj  depends  on  M,  but  the  angular  dependence  on  M  enters  in  higher 
order  terms  of  <|i.  The  center  of  the  elastic  field  is  located  at  the  moving 
point  E. 


Since  the  stresses  o  and  a  are  continuous  at  ij/  *  0,  we  find  from 

x  y 


(3.29a,b)  and  (3.36)-(3.37) 


(2*R  *  “  2R  *  “  k* 


<sr>%  “ 

p  p 


2k 


(3.41a,1 
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where  R  ■  Rp  defines  the  radius  of  curvature  of  the  elastic-plastic 

boundary,  at  least  for  small  values  of  ip.  From  (3.41a,b)  we  obtain 

n/R  ■  2/3.  Substitution  of  this  result  in  (3.41a)  yields 
P 

(2^R~)iS  *4  =  ¥  ’  or  Ro  =  (4/9)  (KI/k)2/2n  (3.42a,b) 

P  P 

For  small  values  of  y/x  and  we  next  consider  the  continuity  of 

t_ ,  and  o„ .  We  use 
R<p  R 

tr^  =  (cos2i/)  -sin2ij/)  o  +  (a  -o^) sini/icos^  (3.43) 


a_  *  2a  sin<J«cosi p  +  o  cos2tfi  +  a  sin2i p  ,  (3.44) 

Rxy  x  y 

in  conjunction  with  (3 . 36) — (3 . 38)  and  the  stresses  in  the  plastic  zone. 

It  may  be  verified  that  is  continuous  to  first  order  in  i (>  .  The  stress 

oR  is  continuous  to  order  unity  by  virtue  of  Eq.(3.42).  Collecting  terms 

to  order  y2  and  \[>z  yields  the  relation 

3k 


x  -  xp:  p2 


4R2 


(3.45) 


Next,  we  consider  the  continuity  of  the  particle  velocity  at  small 
values  of  ty.  We  use 


vco8i l>  -  uslnip 


(3.46) 


u  *  vsintp  +  ucosifj 
R 


(3.47) 
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(3.48) 


Since  [u  ]  *  0  and  [u  ]  *»  0  ,  we  find  from  (3.46)-(3.47) 

R 


(vjcosij;  -  [u]sini//  *  0 
[v]sin^  +  [u]cos<j>  ■  0  (3.49)  I 

Substitution  of  (3.11)— (3.12)  into  (3.48)-(3.49)  and  collecting  terms  of  j 

the  same  order  in  y  yields 

[uQ]  *  0  ,  [Vj]  -  0  [u2]  =  0  (3.50a,b,< 

By  the  use  of  (3.39)-(3.40) ,  we  then  obtain  ; 


x=x  : 
P 


u  ■  a 

2  8ER2 


(3.51a,b) 


where  (2.16)  has  been  used.  For  the  steady-state  problem  (3.51a,b)  imply: 

dv. 

x-x  :  v,  -  -3k/2ER  ,  -  3(2+v)k/4ER2  (3.52a,b) 

P  lx  p  dx  p 

It  appears  to  be  difficult  to  solve  (3.33)  and  (3.34)  rigorously.  Just 

as  for  the  Mode-Ill  case*  an  asymptotic  solution  for  small  values  of  x  can* 

however,  easily  be  obtained  by  considering  solutions  of  the  form 


p2  s  p2/xZ»  vix  s  Vlx^x  (3.53a,b) 

The  appropriate  constants  follow  from  (3.33)-(3.34)  as 

P2  -  -  (3k/2) (1-M) ,  Vlx  =  (-3k/E)(l-M)/M  (3.54a,b) 

The  corresponding  strain  is 


In  the  limit  M  •+■  0,  p2  reduces  to  the  quasi-static  solution  given  by 
Eq. (3.18b),  but  becomes  singular. 

The  similarities  in  the  structure  of  the  equations  for  the  Mode-Ill 
and  Mode-I  plane  stress  cases  suggests  an  approximate  approach  to 

(3.33) — (3.34)  similar  to  the  one  used  for  solving  (2.18)  and  (2.19).  Thus, 
if  an  acceptable  approximation  to  p^  would  be  available  a-priori,  then 

(3.34)  would  be  a  linear  ordinary  differential  equation  for  v^.  A  first 
approximation  to  p^  is  provided  by  the  asymptotic  expression  (3.52a).  This 
expression  has  the  correct  limits  at  M  •  0  (quasi-static  case)  and  M  =  1. 

It  may,  however,  be  expected  that  a  better  approximation  will  be  obtained 
by  adding  a  constant  term  ,  and  use 


P  -  -  |k(l-M)f—  +  — ] 

2  2  Lx2  2xd 


(3.56) 


The  second  term  is  chosen  so  that  (3.33)  is  satisfied  up  to  order  0(M) 


near  x  *  x^.  It  is  noted  that  (3.56)  is  completely  analogous  to  (2.40). 
Substitution  of  (3.56)  into  (3.34)  yields 
d2v„ 


1  - —  -  U-M)  [—  +  —1  v  1 

2  dx2  Lx2  2x2  ■*  lx  Ex' 


(3.57) 


An 


expression  for  x  is  obtained  by  enforcing  the  condition  (3.45)  on  p2: 


[ (1-M) (24M) 


(3.58) 


A  solution  to  (3.57)  is  obtained  by  using  a  perturbation  solution  which 
neglects  terms  of  order  0(M2).  By  integrating  the  result,  the  strain 
e  ■  v.  it  obtained  as 

y  1 
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(3.59a,b) 


c 

y 


^y^PB 


+  eP(x/x  ,M), 

y  P 


where  (e  )D_  =  (k/E)(2-v) 

y  r  D 


3(1-M) 

M 


L  p  \  p 


Also 


a,  = 


i[l  +  (9-8M)1*]  ,  a2  =  |[1  -  (9-8M)*5] 


M(l-M) 


M(l-M) 


2  h 

2[2-(9-8MK] 


0 

C 


A,  = 


2[2+(9-8M)  2] 

.  i]  (.2«262«62,+  P^ossl  .  ]  (1+8 : 

P  P 

-  tFT*1  -  ¥£>  ]  <«1>-  M  «v,  w 


A  *  .(a2-a1)(l+B1+B2+8182)  +  2(8^) 


C 


Here  x  /R  is  given  by  (3.58) 

P  P 

Equation  (3.60)  reduces  to  the  quasi-static  solution,  which  has  been 
given  by  Achenbach  and  Li  (6,Eq.(64),  provided  that  M«.n(x/x^)  <<  1.  In 
the  limit  x  0,  (3.59)  reduces  to  (3.55). 

Numerical  results  for  as  given  by  (3.59)  are  shown  in  Fig.  6,  for 
v  -  0.3  and  M  ■  0.1,  0.3,  0.5.  The  quasi-static  solution  has  also  been 
shown  in  Fig.  6. 


•  60) 

.61) 

.62a,b) 

(3.63) 

(3.64) 

.65) 
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Finally,  just  as  for  the  Mode-Ill  case,  we  apply  the  crack  growth 
criterion  of  critical  plastic  strain  to  determine  the  value  of  that 
would  be  required  for  crack  growth  at  a  given  value  of  M.  For  a 
stationary  crack  the  quasi-static  plastic  strain  follows  from  the  results 
of  Refs. (5)  and  (6)  as 


Vy  -Tc2<f->!  -!<H  • 


(3.66) 


where 


b2  =  fg  Kk+5)  +  (k+1) /2] (E/p) 


(3.67) 


C2  =  16  1"(k+5)  +  (<+1)2/2] (E/p)  (3.68) 

The  constant  k  is  defined  as  k  =  (3-v)/(l+v).  Note  that  properly 
vanishes  at  the  elastic-plastic  boundary. 

Now  suppose  that  the  normalized  critical  strain 

'  <#cr'VPB  »•«) 

is  reached  at  x  =  x^,  for  a  value  of  =  K^.  The  corresponding  value  of 
Xp  is  given  by  Ref . (6,Eq. (57))  as 

yf;  <Kic/k)! 

A  cubic  equation  for  x,/x  follows  from  (3.66).  The  relevant  real- 

r  pc 

valued  root  is 


x,/x  *  S 
f  pc 


(3.71) 


where 
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(3.72) 


s  =  (b2/c2+a)1/3  -  (a-b2/c2)1's 


*  “  1<B2/C2>*  +  [ 


1-f,  +  (E/k)ct/3^h 

c:  J  t 


Next  we  turn  to  (3.60),  and  we  compute  e^x/x^.M)  at  5  =  x^/x^.  Since 
Xp  is  given  by  (3.58)  we  have 


E - fcs  <KIc^> 


(1-M;  (2+M) 

where  (3.42b)  has  also  been  used.  The  crack  growth  criterion  now  yields 


ef  ey^,M)/(£y)pB  * 

where  the  function  form  of  is  given  by  (3.60).  Equation  (3.75)  has 
been  used  to  plot  KjVK^  versus  M  for  three  values  of  in  Fig.  7. 


(3.73) 


(3.74) 


(3.75) 


ft 


Appendix  A 

Conditions  at  a  Fast-Moving  Surface  of  Strong  Discontinuity 

In  a  recent  paper  Drugan  and  Rice  (16  )  have  shown  that  all  stress 
components  are  continuous  across  a  quasi-statically  moving  surface  of 
strong  discontinuity  in  an  elastic-plastic  solid.  They  also  showed  that 
the  only  components  of  strain  which  may  suffer  discontinuities  across 
such  a  surface  are  the  plastic  components  which  do  not  deform  elements 
in  the  plane  of  the  surface,  and  these  strains  may  be  discontinuous  only 
if  the  stress  state  at  the  moving  surface  meets  specific  conditions. 

In  this  Appendix  we  attempt  to  extend  some  of  the  results  of  (16) , 
to  the  case  that  the  surface  moves  so  fast  that  dynamic  effects  must  be 
taken  into  account.  The  results  will  serve  to  establish  the  conditions 
at  the  leading  edge  of  the  elastic-plastic  boundary  of  the  plastic 
loading  zone  ahead  of  a  rapidly  propagating  crack  tip,  particularly  in 
the  immediate  vicinity  of  the  crack  line. 

The  propagating  elastic-plastic  boundary  is  denoted  by  I,  see 
Fig.  1.  The  surface  propagates  with  velocity  V  in  the  direction  of 
the  normal  ^ .  The  coordinate  system  i ®oves  with  the  surface 
l.  A  discontinuity  of  a  field  quantity,  say  g(£^»£2>  Cj**)  is  denoted 
in  the  usual  manner  by 

fig]]  =  g+  -  g  »  (A.l) 

where 

*  =  g(5l,C2,C3,ta  +  &)  (A,2) 
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where  A  >_  0,  and  t  is  the  time  at  which  E  arrives  at  a  particular 
material  point.  In  the  sequel  Latin  indices  i,j,k  have  range  1,2,3. 
Greek  indices  a, 3  have  range  2,3  only  and  thus  refer  to  tensor  com¬ 
ponents  in  planes  parallel  to  planes  that  are  tangential  to  E. 


The  impulse-momentum  relation  yields  across  E: 

llo^n  =  -  pv  [[u  ]]  (A. 3) 

where  p  is  the  mass  density.  By  virtue  of  displacement  continuity  we 
have 

t[uin  -  0  (A. 4) 

It  follows  from  (A. 4)  that,  see  e.g.  Hill  (17), 

[  (3ui/acctl  ]  =  o  (A. 5) 

[[uj]]  «  -  vuauj/a^n  (a. 6) 

By  combining  (A. 3)  and  (A. 6)  there  results 

Cto  II  =  pv^fauj/Hj]  (A. 7) 

Equations  (A. 5)  -  (A. 7)  have  some  implications  for  discontinuities 


in  the  components  of  the  small-strain  tensor,  =  ■j(3u^/8£j  +  au^/Df^). 
We  find 

[[en]l  -  [0^/3^]]  -  l[on]]/pV2  (A. 8) 

[[elaJ3  "7  r [3ua/3ei] ]  -  j[[ola]]/pV2  (A. 9) 
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(A. 10) 
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where  is  understood  to  be  a  stress  state  at  or  below  yield  for  all 
states  along  the  strain  path  from  e+  to  e  .  In  the  inequality  (A. 17), 
WP  is  the  plastic  work  accumulated  discontinuously  at  a  material  point 
due  to  the  passage  of  >. : 


W 


/, 


-P- 

ij 

P+ 

ij 


0  deP 
ij  ij 


(A. 18) 


Subsequent  considerations  are  for  the  special  cases  of  anti-plane 
shear  and  generalized  plane  stress. 

Anti-plane  shear.  This  case  is  defined  by  u^=u  ,t)^  0»  *  u^  =  0. 

The  relevant  relations  between  the  strain  and  stress  increments  across 
I  follow  from  (A. 14)  and  (A. 15)  as 


de?.  -£)do13 


13  2\ 


PV“ 


rfe 


23 


"  2y  dp23 


Substitution  of  (A. 19)  and  (A. 20)  into  (A. 18)  yields 


(A. 19) 


(A. 20) 


”P--  |sl)  +  °13)[lE13n  - 
For  we  now  choose 

o  -  o  + 

°13  “  °13  ’  °23  “  °23  * 

The  inequality  (A. 17)  then  yields 


>  0 


(A. 21) 


(A. 22) 


(A. 23) 
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If  we  restrict  our  attention  to  the  sub-sonic  case,  for  which 


tjV2/v  <  1,  it  is  evident  that  (A. 23)  can  be  satisfied  only  by 

[[o13l]  -  0  and  [[a^]]  =  0  (A.24a,b) 

It  remains  to  verify  that  the  stress  state  (A.22a.b)  is  sub-yield. 

We  should  have  (a^)2  +  (0^3) ^  _S  h2 .  Since  the  yield  condition  is 
satisfied  at  the  -  side  of  I  ,  we  do  have  (c^)2  +  ^°23^2  “  k2 .  ®y 
eliminating  (o^)2,  the  requirement  that  the  stress  state  (A.22a,b)  is 
sub-yield  may  then  be  written  as 

(023)2  "  (°23)2  -  °  (A. 25) 

For  example,  (A. 23)  is  satisfied  if 

°23  >  0  »  °23  >  0  and  °23  -  °23  (A. 26) 

This  is  the  case  that  generally  applies  at  an  elastic-plastic  boundary 
ahead  of  a  crack  tip  propagating  in  Mode-Ill. 

On  the  crack  line  o*3  ■  0^3  =  0,  and,  the  expression  for  Wp  given 
by  (A. 21)  simplifies  to  WP  -  -  (a^  +  a23^e23^'  By  taking 

°23  *  °23  ^  definition  sub-yield),  we  then  obtain  that 

-(l/2y) [ to23l ]z  >  0  ,  which  implies  that  on  the  crack  line 

t[o23l]  -  0  ,  (A. 27) 

without  additional  conditions. 

Generalised  plane  stress.  This  case  is  defined  by 

®33  =  °22^l’^2,t^  ^  ^  0,  ^ 

The  relevant  relations  between  the  strain  and  stress  increments  follow 

from  (A. 13)  and  (A. 14)  as 
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(A. 28) 
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de 


22 


1  .  ,  v  . 

"  E  d°22  +  E  dall 


(A 


dti2  •  T  ^  *°u 


(A 


Substitution  of  (A. 28)  -  (A. 30)  into  (A. 18)  yields 


W* 


-  2(an+0n)  ^ell51  '  2(°22'H,22)[EeP22^ 


(G12'KJi2)[le12n  (A 


For  we  choose 


'11 


all*  °22 


°22  *  °12  °12 


(A 


The  inequality  (A. 17)  then  yields 


If  we  restrict  our  attention  to  p V2/y  <1,  then  (A. 33)  can  be  satisfied 
only  if 

[[au]]  -  0  .  [[a22]]  -  0  and  [[o^l]  -  0  (A 


The  results  (A.34a,b,e)  hold  if  the  stress  state  (A. 32)  is  indeed 
sub-yield,  i.e.,  if 

(a“i>2  +  (<?22^2  ”  ciia22  +  3(o12^Z  -  3k2  (A 

Since  the  stress  state  o^,  o ^  and  o”2  satisfies  the  yield  condition, 

(o^)2  +  3(o”2)2  can  be  exPre88ed  in  terms  of  (c”^)2  and  ail°22* 
Substitution  of  that  result  into  (A. 35),  reduces  that  condition  to 


.29) 

.30) 

.31) 

.32) 

.33) 

.3*a,b 

.35) 


(A. 36) 


+  o 


22 


°ll)[t°22 


]]  <  0 


This  equation  is  satisfied  if  either 


[[o22]]  <_  0  and  ou  <_  o+2  +  °22  (A.37a,b) 

or 

[[o22n  1  0  and  _>  0^2  +  o~2  (A.38a,b) 

Of  interest  in  the  present  paper  are  discontinuities  across  the 
elastic-plastic  boundary  near  the  crack  line,  for  Mode-1  crack  propa¬ 
gation  in  generalized  plane  stress.  For  that  case  we  have  near  the 
crack  line  o+ 2  *  ®  and  022  >  0.  We  also  have  o^  *v<  k  and  022  ^  2k. 

Hence  (A. 3 7b)  is  satisfied.  We  will  generally  also  have  that 
0^2  —  °22  *  hence  f[°22^  i.  0*  and  (A. 37a)  is  satisfied.  Thus,  (A. 32) 
is  an  acceptable  sub-yield  stress  state,  and  the  results  (A.34a,b,c) 
are  valid.  Note  that  on  the  crack  line  -  0*  and  only  [[0^]]  *  0 
and  [[c22]]  *  0  are  relevant. 


32 


Appendix  B 

Exact  Solution  to  Eqs.(2.18)  and  (2.19) 

First  we  introduct  new  variables 

T  =  -Tl/(k/2Rp)  ,  f  =  Yl/(k/pRp),  X  =  x/Rp 
From  Eq. (2.18)  we  then  obtain 

r  -  —  [(i-m2)f  +  ~\ 

M2 

where 

F  =  (l/T2)dT/dX 

Substitution  of  (B.2)  into  the  dimensionless  form  of  (2.19)  yields  by 
the  use  of  (B.3) 

2«-«2>(FitTI)  +2i  +  il2-° 

Equation  (B.4)  can  be  rewritten  as 


2M  — y  = 


/  1-M _ l-m  \ 

'  F  +  F  + 


dF 


2(14M) 

Integration  of  (B.5)  gives 


2 (1-M) 


1-M  1  1+M» 

2M£nT  ■  in  jc|2F+-j^j|  /|2F+Y^jl  j 

where  C  is  an  integration  constant. 

At  X  ■  x  /R  (the  elastic-plastic  boundary)  we  have 
P  P 

T  *  1 


(B.la.b.c) 


(B.2) 


(B.3) 


(B.4) 


(B.5) 


(B.6) 


(B.7) 
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Fig.  3:  Comparison  of  exact  and  approximate  expressions  for 

»  1  +  Low*^  curves  are  approximations 

^  r  y 

according  to  Eq.(2.41). 
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Position  of  elastic-plastic  boundary  versus  M. 
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